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Abstract 



ON ■ We study the cosmological gravitino production in gauge mediated super- 

symmetry breaking models, while properly taking into account the existence 

^^ \ of the messenger mass scale. It is found that for sizable parameter range of 

the model the messenger sector contribution leads to more stringent upper 
bound on the reheat temperature obtained from the condition that the uni- 
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verse should not be overdosed by relic gravitinos. However it turns out that in 
JL , the limit of relatively low messenger scale and large gravitino mass, the relic 

D \ gravitino mass density can be smaller than the critical density independently 



^ 
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of the reheat temperature, suggesting the possibility in this limit to have a 
high reheat temperature without the dilution of gravitinos at late time. 
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I. INTRODUCTION 

Supersymmetry (SUSY) is an elegant solution to the mass hierarchy problem of the 
standard model |Ij|. But it must be spontaneously broken because we do not see degenerate 
pairs of particles and their superpartners. Within the minimal supersymmetric standard 
model (MSSM), it seems unavoidable to introduce explicit supersymmetry breaking through 
soft breaking terms. Such soft breaking terms might be a consequence of spontaneous 
supersymmetry breaking which is transmitted to the MSSM sector by certain mediating 
interactions. In this regards, there are two classes of models which are distinguished by the 
nature of mediating interactions: gravity-mediated models [[IJ and gauge-mediated models 

Cosmological consequences of supersymmetry breaking models have been also discussed 
to some extent. One of the major issues has been the relic abundance of the lightest su- 
persymmetric particle (LSP). In gauge-mediated models, one interesting feature is that the 
gravitino is very light and thus becomes the LSP. For a light gravitino, its interaction is dom- 
inated by the longitudinal goldstino component whose interaction strength is proportional 
to I/7713/2 and thus gets stronger as the gravitino mass 7713/2 becomes smaller ||. Then the 
considerable amount of light gravitinos can be produced from the decay or collisions of the 
MSSM particles, even in the case that the gravitino number density was negligible just after 
the reheating phase of inflation. Not to overdose the universe, this imposes a strong con- 
straint on the reheat temperature depending upon the gravitino mass, unless the gravitinos 
produced in this way were diluted by a late time entropy production ||||. 

In the previous studies j|||, only the gravitino productions by the MSSM particles 
were considered under the assumption that the local goldstino couplings to the MSSM 
fields are valid up to energy scales around the reheat temperature Tr. However it was 
observed by one of the authors (T.L.) that gravitinos decouple from the MSSM fields at 
temperatures above the messinger scale at which the gravitino interaction to the MSSM 
particles is induced radiatively 0. In this paper, we reexamine the cosmological gravitino 
production while properly taking into account the existence of the new mass scale in gauge- 
mediated models, i.e. the messenger scale Mx, which is deeply involved in the problem. 
We first note that the goldstino production due to the MSSM particles at high temperature 
T > Mx is suppressed by M x /T 2 when compared to the naive extrapolation of the low 
temperature results at T < Mx- However, in this high temperature range the messenger 
sector contribution is much larger than that of the MSSM particles because the messenger 
particles have stronger couplings to the goldstino. As a result, for much of the parameter 
range of (M x , 7n 3 / 2 ), the messenger sector contribution leads to more stringent upper bound 
on the reheat temperature which is obtained from the condition that the universe should 
not be overdosed by relic gravitinos. The results are depicted in Fig. 2 for several values of 
the messenger and MSSM gaugino masses. 

An interesting feature of the messenger sector contribution is that the relic gravitino 
number density to the entropy number density ratio (I3/2 = n 3/2(T)/s(T)) behaves such 
as dY 3 / 2 oc dT / {rn^^T 71 ) (n > 2) for M x < T < T R , indicating that the present value of 
Y3/2 is mostly induced at T ~ Mx, not at the highest available temperature T ~ T R . As 
a result, in gauge-mediated models with Mx ~ T R , the gravitino abundance is determined 
by the values of m 3 / 2 and M x , and is independent of T R . This means that in the limit 



of relatively low messenger scale and large gravitino mass, the relic gravitino mass density 
can be smaller than the critical density independently oJTr. Of course, this would be true 
only when other more model-dependent contributions, e.g. those from the supersymmetry 
breaking sector, are negligible, which would be the case if the typical mass scale of the 
supersymmetry breaking sector is bigger than Tr. The results of Eq.([L9[) and Fig. 2 show 
that for the messenger mass Mx = 10 5 ~ 10 6 GeV and the gravitino mass 7713/2 ^ 100 
MeV, we can have Tr much higher than the previously obtained bound without assuming 
the dilution of gravitinos at late time. 

II. GOLDSTINO COUPLINGS IN GAUGE-MEDIATED MODELS 

Gauge mediated supersymmetry breaking (GMSB) models consist of the observable 
(MSSM) sector, the messenger sector, and also the supersymmetry breaking sector 0. Su- 
persymmetry breaking sector breaks supersymmetry dynamically at certain high energy 
scales and delivers supersymmetry breaking effects to the messenger sector through mes- 
senger gauge interactions or other interactions [§J§. Supersymmetry breaking in the 
messenger sector is finally transmitted to the MSSM sector via the standard model gauge 
interactions and generates the soft supersymmetry breaking masses of order 10 2 ~ 10 3 GeV. 
After integrating out the supersymmetry breaking sector, the effective superpotential is 
given by 

W eS = \X<S><S> c + AW(X,Y i ), (1) 

where X and Yi are the standard model singlets and the nonperturbative effects of the 
supersymmetry breaking sector are encoded in AW. Here we consider the simplest mes- 
senger sector, consisting of single flavor of chiral superfield $ + $ c transforming as 5 + 5 of 
SU(5i). By minimizing the effective potential from W e g, we can determine the vacuum ex- 
pectation values (VEVs) (X), (Y*) and also (F x ) = (dW eS /dX), (F*) = (dW^/dY^ which 
measure the size of supersymmetry breaking. The VEV (X) gives a supersymmetric mass 
X(X) = Mx to the messenger $ + $ c , while (F x ) leads to the supersymmetry breaking mass 
splitting in the scalar masses of <J> + $ c . Supersymmetry breaking in the messenger sector is 
transmitted to the MSSM particles via the standard model gauge interactions at one-loop 
or two-loop levels, yielding the supersymmetry breaking mass splitting in MSSM particles: 

m A ~ (a/Air)F x /(X), m) ~ («/4vr) 2 (F x /(X)) 2 , (2) 

where irt\ and m? denote the soft gaugino mass and scalar mass, respectively, and a is the 
MSSM fine structure constant. 

The spontaneous breakdown of supersymmetry leads to a massless spin-1/2 fermions, 
the goldstino. The goldstino field is a linear combination of the fermion components of X 
and Yi weighted by F x and F iy i.e., 

X = b > ( 3 ) 



where F = J \Fx\ 2 + Y^i\Fi\ 2 ■ When global supersymmetry is promoted to local super- 
symmetry (supergravity), the goldstino is absorbed into the gravitino as the longitudinal 
components and the gravitino acquires a mass 



'" S/2 = Tufr (4) 

with the vanishing cosmological cons tant assu med. 

At energy scales below y/F ~ Jm^^Mp but above the messenger mass scale Mx = 
X(X), the goldstino x directly couples to the messenger fields, which can be read off from 
the superpotential ([!]) using the relation ipi = (Fj/F)x H , 
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where <fr + <P C and ip + ip c denote the scalar and fermion components of the messenger super- 
fields $ + $ c . At these energy scales, there is no direct local couplings of the goldstino to 
the MSSM particles at lagrangian level. However at lower energy scales below Mx, one can 
show H that integrating out the messenger particles through one and two loop diagrams 
leads to the following effective lagrangian P, |lO|jrT| : 



2 2 

Cmssm = ^ mf fxf* + ^ W^XF^ - r - 9mx /*/xA + h.c, (6) 

y3m 3 / 2 Mp AyQm^^Mp yQm^^Mp 

where (/, /) and (A, F^ v ) denote the MSSM matter and gauge multiplets, respectively. We 
stress that the above local goldstino couplings to the MSSM fields are valid only at energy 
scales below Mx- When computing the amplitudes with an external energy E > Mx (but 
still below y/F), one has to use (pj), not @, even when the external particles involve only 
the goldstino and MSSM particles. Note that the effective lagrangian @ summarizes the 
leading order results of the low energy amplitudes expanded in powers of E/Mx- 

For the soft masses (|2|) presumed to be of order 10 2 ~ 10 3 GeV, goldstino couplings are 
described essentially by the two unknown model-dependent parameters, m 3/ / 2 and Mx- In 
the next section, we will compute the relic gravitino mass density as a function of 777.3/2, 
Mx, and also the reheat temperature Tr. The size of 7713/2 in GMSB models has a wide 
range of model-dependence, which can range from 0(1) eV even to 0(10) GeV. Moreover, 
depending on how the messenger scale is generated and also on how supersymmetry breaking 
is transmitted from the supersymmetry breaking se ctor to th e messenger sector, it is possible 
to have M x and \/Fx much smaller than y/F ~ Jm^^Mp. (Note that F x /M x is fixed to 
be of order Aitm\/a ~ 5 x 10 4 GeV.) This is an interesting possibility in connection with 
allowing a large reheat temperature, which will be discussed in the next section. 



III. COSMOLOGICAL GRAVITINO PRODUCTION 

At first, let us estimate the temperature range for which the gravitinos are thermalized. 
For the gravitinos to be thermalized, their reaction rate is required to be larger than the 
expansion rate of the universe. The expansion rate of the universe is measured by the 
Hubble parameter H(T) = (7r 2 g*/90) 1//2 T 2 /Mp for the radiation dominated universe. For 



the temperature below the supersymmetry breaking scale yF ~ Jm 3 / 2 Mp, the reaction 
rate of the gravitinos is roughly given by 

r(T) ~ ZMmiT), (7) 

i 

where (07) is the thermal averaged cross section producing the gravitino and ni(T) is the 
number density. For T < Mx, the number density of the messenger particles is Boltzmann- 
suppressed and then we have to consider only the contributions from the MSSM particles, 
while for T > Mx, the contributions from the messenger particles are dominant. The 
resulting reaction rate is given by 



r(T) 
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where £2 and £4 are constants of order unity. (For £ 2 , £4 and the cross sections, see the 
discussions below.) Then from T(T) > H(T), we obtain the following temperature range 
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for which the gravitinos are thermalized. Once the gravitinos were thermalized, the gravitino 
mass density to the critical density ratio at present is given by [|T2 



a^ ~ (if) • do) 

and so 7713/2 should be less than about 1 keV not to overdose the universe now. (Here 
h = ifo/(100km/sec/Mpc) for the present Hubble expansion rate Hq.) 

The simplest way to avoid the result ([H]) is to assume that gravitinos had never reached 
to thermal equilibrium after the inflation, which would allow m 3 / 2 ^ 1 keV. Even in this 
case, significant amount of gravitinos can be produced by the decays or collisions that occur 
in the thermal bath of the early universe subsequent to inflation |TBj , which usually leads to 
an upper bound on the inflation reheat temperature. We can calculate the number density 
of such gravitinos using the Boltzmann equation 

h 3/2 + 3Hn 3/2 = J2( r (i->x+-)) n i + J2( a (.i+J^x+-) v ) n i n J H • ( n ) 

/ i,j 

Introducing Y 3 / 2 {T) = n 3 / 2 {T) / s{T) where s(T) = (27r 2 /45)^ s T 3 is the entropy density, 
and changing the variable from the time to the temperature, we can integrate the above 
Boltzmann equation to obtain 



fn 3 /2\ = [ Tr , r E/(rV^ x+ ...))n J + Ei,j(v(i+J^ x +-)v)ninj + 
V Jo k s{T)H{T)T 



Y — ( 3/2 ) = dT (±^x±zLL_L 21 ^i,J\ u ^+J^x+-) L, / l ^"'J 2 (12) 



Here Tr is the reheat temperature after the inflation and the RD era formula —dt/dT = 
1/HT is used since most contributions come during the radiation dominated era. The 
evaluation of the integral of ([12]) for given r(/_ +x+ ...) and cr(/ + j_ >x+ ...) will be explained in the 
appendix. The gravitino mass density produced by the decays or collisions in the thermal 
bath is then given by 

n 3/2 h 2 ~ 2.8 x 10 8 F3/2 (^) • (13) 

The processes which dominantly produce the gravitinos depend on the temperature of the 
thermal bath. For T < Mx, MSSM particles are the dominant source of light gravitinos. 
From the effective lagrangian (§), T^ x ^ y+X) and a( x+y ^ z+x ) where x, y, z are the MSSM 
particles can be calculated. When summed over all possible channels, the results are M 

Y, T (^v+x) = £1 —2 A A/f2 > ( 14 ) 

x,y rn 3/2 1V1 P 

7; 2J a {x+y-z+ x ) = 6 2 /if 2' ( 15 ) 

z i,3/,z rn 3/2 1V1 P 

where £1 and £ 2 are of order 1. We then find for m x < T R < Mj 

"^ ^ (l^vf H (lo^v) 3 + 2 ^ (lo^v) 2 (lolv)] < 16 > 

For T > Mx, the effective lagrangian (|]) which gives the cross section (fTfj|) is not valid 
any more. When computed using the correct lagrangian (|5]), the cross section a^ x+y ^ z+x ){s) 
falls like 1/s for s/s > Mx- As a result, the gravitino production due to the MSSM particles 
at temperature T > M x is suppressed by M x /T 2 when compared to the naive extrapolation 
of the result obtained from the effective lagrangian (^). Now the gravitinos are produced 
dominantly from the processes involving the external messenger particles which couple to 
the goldstino more strongly. The relevant decay widths and scattering cross sections are 
given by 

&*-**'«> -t*\^) 3^MT 

/47TX 2 m\M\ 
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where X, X' represent the messenger particles transforming as 5 + 5 of £77(5) and 



(17) 



6 = 7~> ^ 4 = «i + 3a 2 + 8a 3 . (18) 

It turns out that the most important contribution comes from the messenger particles at 
T ~ Mx, and thus we need a more careful evaluation of the production cross sections 
for s = 0(M X ). Now the relic gravitino abundance from the messenger particles can be 
obtained using the equations in the appendix. For T R > M x , we find 
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• (19) 



In Fig. 1, we plotted VL 3 / 2 h 2 including both the MSSM contributions and the messenger 
contributions as a function of the reheat temperature Tr for typical values of the gravitino, 
messenger and MSSM gaugino masses. It shows that f2 3 / 2 for Tr < O.lMx is from the 
MSSM contribution, while 1)3/2 for Tr > O.lMx is essentially from the messenger sector 
contribution. Note that f2 3 / 2 is almost independent of Tr for Tr > Mx as summarized in 
(|19|) . This is because Y 3 / 2 = n 3 / 2 {T) / s{T) produced by the messenger particles behaves such 
as 



dT 

J 3/2 J 



dY3/2(X m^ (n " 2) ' (20) 



and so the present value of Y 3 / 2 is mostly induced at T ~ Mx, not at the highest available 
temperature Tr. This feature requires a more careful treatment of the integral (|T2| ) whose 
integrand takes a nontrivial form for T ~ Mx- 

Now we can obtain an upper bound on Tr as a function of the gravitino mass by requiring 
^3/2^ 2 < 1- We present the results in Fig. 2 for several values of the messenger and MSSM 
gaugino masses. Summarizing the results, for the gravitino mass range 
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we find the practically same result as the previous one obtained considering just the MSSM 
contributions |4j]. For larger m 3 / 2 , the messenger sector contribution leads to a significantly 
stronger bound on the reheat temperature: 



T R < O.lMx (22) 



for the gravitino mass range 



0.1 ( ^) k^y MeV < m 3/2 < 10= (^) f ^ MeV. (23) 



10 3 GeV/ V10 5 GeV/ " " /z V10 3 GeV/ Vl0 5 GeV7 

However no constraint is imposed on Tr for the gravitino mass range 

\ 2 f M x \ 



m 3 /2 



fc'GravMu&O™- < 24 > 



suggesting that in gauge mediated models with a small messenger mass (Mx = 10 5 ~ 10 6 
GeV) and also a relatively large gravitino mass (7713/2 ^ 100 MeV), it may be possible to 



have a rather high reheat temperature without the late dilution of relic gravitinos. Of course, 
this is possible only when other more model-dependent contributions, e.g. those from the 
supersymmetry breaking sector, are small enough, which would be the case if the typical 
mass scale of the supersymmetry breaking sector is bigger than T R . 
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APPENDIX 



In this appendix, we present the formulae for the integral of (0) when the decay width 
T(E) and the cross section a(s) are given. The decay and scattering terms in the Boltzmann 
equation can be written as 
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where V is the decay width in the rest frame and K% is the Bessel function. Inserting these 
into flT2| ) and changing the order of integration, we obtain 
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where g = 135v / 10/(27r 3 5fJ /2 ) and t R = m/T R in Y Acc . Ay and (mi + m 2 )/T R in F S cat- For y scat , 
the Maxwell-Boltzmann distribution is used for the simplicity of calculation, which gives the 
correct result up to a factor of order 1 . 

For T R > m, Fdccay is well approximated by 
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For Kjcat, we present two typical examples. First, when mi,m2 -C T R and <x(s) 
constant, Y SCSut is given by 
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For the second example with mi = M x , m 2 = and a(s) of (|I7D , we have 
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FIG. 1. J7 3 / 2 /i 2 as a function of the reheat temperature Tr for the gravitino mass m 3 / 2 = 10 
MeV, the MSSM gaugino mass m\ = 300 GeV and the messenger mass Mx = 10 5 GeV (solid line), 
10 6 GeV (dashed line), 10 7 GeV (dotted line). tt 3 / 2 h 2 is scaled by (10MeV/m 3/2 ) for different 
values of m 3 / 2 ■ 
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FIG. 2. The upper bound on the reheat temperature Tr from fi 3 / 2 h 2 < 1 as a function of the 
gravitino mass 1713/2 in the gauge mediated models. We used the MSSM gaugino mass m\ = 300 
GeV in (a) and m\ = ITeV in (b), and the messenger mass Mx = 10 5 GeV (solid line), Mx = 10 6 
GeV (dashed line), Mx = 10 7 GeV (dotted line), respectively. 
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